Damped bounces of an isolated perfect quantum gas 
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The issue of the thermahzation of an isolated quantum system is addressed by considering a 
perfect gas confined by gravity and initiaUy trapped above a certain height. As we are interested 
in the behavior of truly isolated systems, we assume the gas is in a pure state of macroscopically 
well-defined energy. We show that, in general, for single-particle distributions, such a state is strictly 
equivalent to the microcanonical mixed state at the same energy. We derive an expression for the 
time-dependent gas density which depends on the initial gas state only via a few thermodynamic 
parameters. Though we consider non-interacting particles, the density relaxes into an asymptotic 
profile, but which is not the thermal equilibrium one determined by the gas energy and particle 
number. 

PACS numbers: 03.65.-w, 05.30.-d, 05.70.Ln, 05.60.Gg 
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I. INTRODUCTION 

The dynamics of a quantum system under the influence 
of a heat bath has been widely studied in the last decades. 
Under these canonical conditions, the system relaxes into 
thermal equilibrium at the bath temperature, provided 
the coupling to the bath is weak enough [1[. Relaxation 
under microcanonical conditions, assumed in thermody- 
namics, is far less understood. A possible definition of 
microcanical conditions might be an energy-conserving 
coupling to environmental degrees of freedom @, S] ■ ^s 
is well known, the measurement process is convincingly 
accounted for by such a coupling of the measurement ap- 
paratus to its environment: the coherences between the 
measured system eigenstates are generically destroyed 
whereas the populations remain constant [J, la, |0|. It 
therefore seems difficult to regard the evolution under 
such conditions as a relaxation process. However, one 
can wonder whether physically relevant system degrees 
of freedom relax towards their thermodynamic equilib- 
rium values. This question is also of interest for truly 
isolated systems, i.e., in absence of any coupling to envi- 
ronmental degrees of freedom. 

This issue is especially relevant for the n-particle re- 
duced density matrices of a many-body system. Re- 
cently, different boson and fermion systems have been 
studied numerically [7,, iS, [9,, illi] ■ In these works, the sys- 
tem is initially prepared in the ground state of a Hamil- 
tonian different from that governing the time evolution 
and single-particle distributions are evaluated. The dis- 
tributions considered were found to relax into thermal or 
non-thermal distributions irrespective of the integrable 
or non-integrable nature of the system. In [ll| , the Joule 
expansion of an isolated perfect quantum gas has been 
obtained analytically. In this study, the gas is initially 
trapped in a subregion of the entire accessible volume but 
is not assumed to be in the corresponding ground state. 
Following 0, y, [13, [S [3 ) a pure state of macroscop- 
ically well-defined energy has been considered and the 
thermodynamic limit expression for the time-dependent 



particle number density has been derived. Interestingly, 
for this single-particle distribution, the initial gas state 
is fully characterized by a few thermodynamic param- 
eters, such as the gas energy and particle number or, 
equivalently, the corresponding microcanonical tempera- 
ture and chemical potential defined from derivatives of 
the gas microcanonical entropy. 

In this paper, we consider an isolated quantum gas 
confined by a uniform gravitational field and a perfectly 
reflecting wall. The gas is at first trapped above a cer- 
tain height and then left free to fall. Such experiments 
can be realized with laser-cooled atoms and bounces of 
both Bose-Einstein condensates and thermal clouds have 
been observed [1^, [la] ■ We study the case of a perfect 
gas. The model we examine is presented in the next sec- 
tion. In Sec. Illli we discuss the states of macroscopically 
well-defined energy. We show that, for single-particle 
distributions, such a pure state is strictly equivalent to 
the microcanonical mixed state at the same energy. Our 
derivation is applicable to any many-body system. For a 
perfect gas confined by gravity, we are able to derive the 
time-dependent density profile in the limit of a large par- 
ticle number. This is done in Sec. lIVI We obtain damped 
oscillations of the gas density. However, the asymptotic 
particle number density of the gas is far from the thermal 
equilibrium one determined by the gas energy and par- 
ticle number. We then compare the relaxation behavior 
of a truly isolated system to that of a system under the 
influence of an environment. Finally, in the last section, 
we summarize our results and draw conclusions. 



II. MODEL 

We consider a one-dimensional perfect gas confined by 
a homogeneous gravitational field and an infinite poten- 
tial wall. The N particles of mass m constituting the gas 
are described by the Hamiltonian 



Hn = 



dzil}\z) 



-i^'*(=' 



- rngzipiyz) 



(1) 



where ip^{z) creates a particle at position z and g is the 
acceleration of the gravitational field. Throughout this 
paper, we use units in which h — ks = 1- The non- 
interacting Hamiltonian ([T]) can be easily generalized to 
higher dimensions. In the following, we treat in detail 
only the one-dimensional case but we also mention re- 
sults for the 3D case. The Hamiltonian ([T]) can be char- 
acterised by the microscopic energy eg = {mg'^ /2)^'^ and 
length zq = {2m?g)~^/^. For Rb atoms in the Earth's 
gravitational field, eg and zq are, respectively, of the order 
of 10 nK and 0.1 /J,m. The single-particle eigenenergies 
and normalised eigenfunctions of Hq are 



eolflql 



4>q{z) 



A{z/zo + aq) 
'z^A'{aq) 



Q{z) (2) 



where 8 is the Heaviside step function, A{x) — 
J^ dtcos{t^/3 + xt)/TT the Airy function QJi], A' its 
derivative and Ug its zeros. We label these zeros by pos- 
itive integers in the order of increasing absolute values. 
Initially, the gas is trapped above a certain height Z and 
hence described by the Hamiltonian Hz given by the ex- 
pression ll]) with the lower limit replaced by Z. The 
single-particle eigenenergies and eigenfunctions of Hz are 
simply ez,k = mgZ + ek and (/)z,fc = Mz ~ Z). 



III. MICROCANONICAL EVOLUTION 

As mentioned in the introduction, the gas is assumed 
to be initially in a pure state of macroscopically well- 
defined energy E 



i*>= Y. *"i") 



(3) 



\a)eHE 



where \a) refers to the A^-particle eigenstates of the 
Hamiltonian Hz- The Hilbert space Tig is spanned by 
the states \a) corresponding to eigenenergies in the in- 
terval [E, E -f 5E] where 5E is much smaller than E but 
much larger than the maximum level spacing of Hz in 
this interval. The subsequent time evolution of the gas is 
governed by the Hamiltonian H^ which is different from 
Hz- In this section, we first show that, in the large N 
limit, the time-dependent gas density is the same for al- 
most all states ^ and equal to a microcanonical average 
at energy E of the particle number density operator. We 
stress that this result does not rely on the specific form 
of the Hamiltonians iJo and Hz- Our derivation does not 
even have recourse to their non-interacting nature. We 
then discuss the particular Hamiltonian ([T]). 



A. Microcanonical typicality 

The gas density is given by 

p(z,t) = {^e'"«'^\z)^{z)e~'"°'\^)- (4) 



To prove that p is the same for almost all states ([3]), we 
use, following Refs. [lil[ij, [3|i the uniform measure on 
the unit sphere in He 



Ai({*J) 



[D^l) 



■5 1 



\a)enE 



(5) 



where D is the dimension of He- We assume that, for 
large N , this dimension is practically proportionnal to 
5E and that a density of states can thus be defined as 
n{E, N) = D/5E ^. This density satisfies Boltzmann's 
relation 



InD = In [n{E, N)6E] = S{E, N) 



(6) 



where S is the gas entropy. 

We now evaluate the Hilbert space average and vari- 
ance of the number 

l-z+S 

Ns{z,t)= / dz'p{z',t) (7) 

-J z 

following from the normalised distribution ([5]). The 
above length S > can be as small as we want. 
For N ^ 1, the reduced distribution of a fi- 
nite number j of components ^,^(1), . . . , ^^(j), equals 

'n.i=iDe-Kp{-D\^aii)\'^)/'^- With this Gaussian distri- 
bution, we obtain the Hilbert space variance 

i-z+S 



Ns'~Ns = 



1 



E 



\a),\l3)eHE 



dz'p{z',t)\(3) (8) 



where p = exp{iHot)tp''tpexp{—iHot) and ... denotes 
the average with respect to the measure ([5]). Upper 
bounding the sum over the states \f3) G He by a sum 
over all the A^-particle eigenstates |/3) gives 



Ns^ - iv/ < 



z+S 



z+5 



dz' / dz"{p{z',t)p{z",t)) (9) 



where (...) ~ X^Iq^p-h (c^I • ■ • l*^)/^ is a microcanonical 
average at energy E. We reiterate that \a) are the eigen- 
states of the Hamiltonian Hz- For the Hilbert space av- 
erage of the density p, we find (p) . As the density-density 
correlation function on the right side of ([9]) is positive, 
the Hilbert space variance of Ns is lower than N"^ /D. 
Consequently, as D — exp(iS') '^ exp(iV), this variance 
vanishes in the limit iV ^ 1 and hence p{z,t) = {p{z,t)) 
for almost all states ([3]). 

The above derivation is applicable to any single- 
particle distribution (*|e*^"*ctcpe~'^''*|*) where c], = 
J dzxp{z)ip'' {z) and the wavefunctions Xp form a basis 
of the single-particle Hilbert space. Moreover, since the 
only property of the observable A — J ' dz' p used in 
this calculation is (A^) <C exp(A^), the equivalence be- 
tween the average (\I'| . . . |^) and the microcanonical av- 
erage (...) shoud be valid for a large class of few-body 
observables. We remark that this equivalence is obtained 
here for a typical state |^) not for an eigenstate |a). Con- 
trary to Refs. [i3i[3i we have not used any remarkable 
property of these eigenstates in our derivation. 



B. Perfect gas confined by gravity 



For the non- interacting Hamiltonians Hq and Hz, 
it is useful to define the creation operators cj, = 

J dz(l)p(z)ip^ (z) and c^ given by similar expressions with 
4>z,p in place of (j)p. With these definitions, the particle 
number density reads 



p(z,t)=^(4c,)0p(z)</>,(z)e^*(^-- 



where 



E 



p.q 



(10) 



a .. 



\a)/D as shown above. 



Here, an eigenstate |a) corresponds to a set of occupation 
numbers {uk} obeying J2k "*: = ^- ^'^^ fermions, Uk is 
restricted to the values and f . The Hilbert space He 
is spanned by the states \a) satisfying E < J^k "^kf-z^k < 
E + SE. The expression (fTU|) can be further simplified 
using the following standard arguments [l^. The mi- 
crocanonical probability distribution of the occupation 
number Uk is P^nk) = n{E - nkez,k,N - nk)/n{E,N) 
where n is the density of states of the Hamiltonian 
Hz — f'Z,kCz k^^z k- This density obeys Boltzmann's re- 
lation ^ with the corresponding entropy S. Expand- 
ing this entropy in the energy nkez,k ^ E and number 
nfc <C iV and taking into account that S* = 5 in the large 
N limit, result in P{nk) (x exp[—nk{ez,k — ^J)/T] where 
the temperature T and chemical potential /i are deter- 
mined by 

^^dES{E,N) , ^^-dNS{E,N). (11) 

The microcanonical averages in (jlO[) can thus be rewrit- 



ten as 



(4^) 



{p\k){q\k) 



Z^e^./T/^zpl 

fc>0 ' ^ ^ 



(12) 



where the upper sign is for bosons and the lower sign is 
for fermions, C — exp[(/i — mgZ)/T] and 

A{Z/zo + aq) 



{<l\k) 



dz<j)q{z)<j)k{z - Z) = 



A'{aq)[aq - ak + Z/zq] 

(13) 
The second equality is derived in the Appendix. 

Other choices than ^ are possible to describe a gas 
state of macroscopic energy E. For example, due to the 
exponential A^-dependence ^ of the density n{E,N), 
one obtains the same gas density p for typical states of 
He and of the Hilbert space spanned by the eigenstates 
\{nk}) satisfying J2k^k = N and J2k'^k<^z,k < E. We 
also remark that, in three dimensions, for a free horizon- 
tal propagation, the number of particles p(z, t)dz between 
z and z + dz is given by pO| and (fT2l) with the thermal 
occupation function replaced by a sum over horizontal 
modes determined by the initial confining potential. If, 
for example, the gas is initially trapped in a square of side 
length y/S ^ (mT)"-*^/^, the thermal occupation factor 
in ini) is replaced by TTOS'rin[l=FCexp(-efc/T)]/27r. In 
this case, the behavior of p is qualitatively similar to that 
discussed in the following. 



< 
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z-Z [ units of T/mg ] 

FIG. 1: Initial particle number density as a function of height 
for fermions, fugacity (" = 10 and temperatures T = 2eo and 
5eo. The dashed line is the large A^ limit. 



IV. DAMPED BOUNCES 

In this section, we derive the large A^ limit expression 
for the time-dependent gas density given by (jlOp and 
P^ . We find that it relaxes to an asymptotic density 
profile Poo(-z)- However, this density poo is not the ther- 
mal equibrium one determined by E and A^. It depends 
not only on these thermodynamic parameters but also on 
the height Z. 



A. Temperature, chemical potential and density 
profile of the initial state 

The temperature T and chemical potential /x are de- 
termined by the gas energy E and particle number A^ via 
pip or via the more convenient relations 



A^ 



fc>0 



/C l«fel 



E 



I/C (l' 



eoy^|a/c|/c flflfc' 



k>Q 

(14) 
where f(^{x) = [exp(a;)/C T 1] ^- We remark that the fu- 
gacity C defined above and the temperature T depend 
only on A^ and E and not on the height Z. As the 
large zeros of the Airy function are essentially given by 
\ak\ — (377/0/2)^/"^, the density of levels e^ = eo\ak\ in- 
creases with k and the above sums are dominated by the 
contributions of zeros |afe| '--^ T/eo which are very close to 
their large-fc estimates. Moreover, for T :s> Cq, the sums 
(fT4| can be approximated by integrals [20| and hence 



A^ = 



3/2 



hMO 



eo 



5/2 



^3/2(0 (15) 



where lu{0 = Jo dxx^ f q{x) / -k . For a Bose gas, these 
expressions hold only for temperatures T above the con- 
densation temperature Tb = eo[Af//i/2(l)]^^^- Below 
this temperature, the right sides of the equalities HI 



with C = 1 give the non-condensed fraction of the initial 
state and the corresponding energy. 

The initial particle number density p{z, 0) = 
J2k>o fci^k/T)(f)k{z — Z)"^ can be evaluated in a simi- 
lar manner as follows. As the Airy function A{x) van- 
ishes rapidly with increasing x, the sum over k is dom- 
inated by the terms \ak\ ^{z — Z)/ z^. For these terms, 
using the large negative x approximation of A{x) and 



(j)k{z-Zf 



i/2|afc|i/4, wefind 

l + sin[4[|afc|-(z-Z)/zo]3/V3] 



(16) 



The contribution of the sine term to (o(z,0) is negligible 
as it oscillates strongly with k and we finally obtain the 
large TV expression 



p(z,0) 



1 



V^K 



dx 



k 



+ 'I^{z-Z)\Q(z~Z) 



(17) 

where A = (27r/mT)i/2 ex {eo/Tf/'^T /mg is the de 
Broglie thermal wavelength. The number of particles be- 
tween z > Z and z -I- d^ is equal to that of free particles 
in a box of volume dz at temperature T and chemical 
potential /i — mgz. This expression does not describe the 
increase of p in the small region z — Z < A but is then 
very accurate, see FigHJ The results shown in FiglT] are 
obtained by numerical summation over the exact eigen- 
states ^ . The characteristic length scale of the particle 
number density is thus essentially T /mg for C < 1- This 
length is of the order of 1 mm for Rb atoms at 0.1 mK 
in the Earth's gravitational field. As C, and T are inde- 
pendent of Z , the thermal equilibrium density profile of 
the gas described by the Hamiltonian iJo, energy E and 
particle number N is given by (J17p with Z — 0. 

In the classical Maxwell-Boltzmann limit, A ^ 
T/mgN, the expression ([T7]) simplifies to p{z,0) ex 
exp[mg{Z — z)/T]. This exponential decrease is always 
valid for large z. For a Fermi gas, in the opposite limit, 
the expressions ([T5|) and P?|) hold for a Fermi temper- 
ature Tp = eo(37rAf/2)2/3 > g^. In this case, ^ be- 
comes p{z,0) ex [1 -I- mg{Z — z)/Tf]^^^. For a Bose 
gas below the temperature Tb, the density profile of the 
non-condensed gas is given by pT|) with C = 1 whereas 
the particle number density of the condensate is propor- 
tionnal to A{z/zo + ai)^ which practically vanishes for 
z > zo = co/mg. For T <Tb, the ratio of the sizes of the 
condensate and thermal cloud is of the order of A^^^/s 
We remark that though (fT7)l with C = 1 diverges in the 
limit mg{z — Z)/T — > 0, the number of non-condensed 
bosons between Z and Z + zq is negligible. 



B. Density relaxation 



The time evolution of the density pO|) can be charac- 
terized by the motion of the center of mass 



zcit) =N^^ / dzzp{z,t). 
Jo 



(18) 
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FIG. 2: Gas center of mass as a function of time for fermions, 
fugacity (^ = 10, temperature T — 5eo and Z = T/mg. The 
dashed line is the large A'^ limit. This approximation is excel- 
lent already for T ~ lOeo- 



As detailed in the Appendix, by evaluating the second 
order derivative of zq and using the property dz(f)q{0'^) = 
{2m'^gy/^, we find 

-^J{Lj)[l~COs{Ljt)] 

J{lo) = A^-i^(ctc,)<5(u;-6p + e,). (19) 

To determine the large N limit of the function J, we first 
rewrite it as 



OO /"OO 



J{uj) = N-^ / dz j dz'n{z,z') 
Jo Jo 



(20) 



x^d{uj-ep + eg)(j>p{z + Z)4)q{z' + Z) 



vi^i 



niz,z') - ^/c(efe/r)0fe(z)0fe(z') 



(21) 



fe>o 



A / '^^•^'^ 



2 , 1^9 

TTX H Z 

T 



„i27r£c(z-z')/A 



The approximate expression for f2 holds for \z — z'| < A. 
For larger |z — z'|, Q, practically vanishes. The func- 
tion ri is simply related to the single-particle density 
matrix at initial time ('0 (z)V'(z')) = 0,{z — Z^z' — Z), 
see ([T71) . Due to the factor ^{z,z'), the main con- 
tribution to the integral (|20p comes from the region 
z ~ T/mg and \z — z'|<A. The sum in the expres- 
sion pop is thus dominated by the terms for which 
ep,eq ~ T and \ap — a^l < (eo/T)^/^^ Pqj- these terms, 
ep — eq ~ 7reo(p — g)(eo/eq)^'2, \Ye use this approximation 
and rewrite the sum in ([^D]) as a sum over q and s — p~q. 
For T ^ Co, the sum over q is well approximated by an 
integral and the second factor in the integrand in (PO)) is 
found to vanish for |z — z'| ^ A. Consequently, the lower 
limit of the integral over z' can be pushed to — cx) and we 
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FIG. 3: Gas center of mass as a function of time in the large 
N limit for fermions and {(,mgZ/T) = (0.1, 1), (100, 1) and 
(10,5). 



obtain 

ZG{t) 



zg{0) 






U-^r^^r'^ (22) 



sin ^/(sTrp 



3Z[(s7r/x)2-l]/^-p I 



1 — cos X\ 1 

2Z 



where Z — mgZ/T. We remark that, sinrilarly to 
Ref. [ll|, it remains a sum over a discrete index in this 
large N expression. For Z ^ T/mg, the time {Z/gy^^ 
is of the ordrer of 10 ms for Rb atoms at 0.1 mK in 
the Earth's gravitational field. To better apprehend the 
characteristics of the many-body case studied here, it is 
instructive to compare with the case of a single parti- 
cle. The average position of a particle can be written 
in the form ((I9l)-(ll0|) with Z = and n{z,z') replaced 
by (p{z)*(p{z') where ip is the initial wave function of the 
particle. This density matrix is very different from Q. Its 
extents along and perpendicular to the diagonal z' ^ z 
are the same whereas fl is almost diagonal which plays 
an essential role in the above derivation. 

In the limit t ^ {Z/gY^^, the contribution of the co- 
sine term in (|22|) vanishes. This expression thus describes 
the relaxation of the center of mass zq to an asymptotic 
value zg{oo). From pT]) and ([M]) . it can be shown that 
zq{oo) — zg{0) — Z/3. A gravitational potential energy 
of only NmgZ/S is converted into kinetic energy whereas 
the potential energy difference between the initial density 
(|17p and the thermal equilibrium one given by ^17\i with 
Z = is NmgZ. The function (j22p is non-nionotonic, 
the motion of the gas center of mass consists of damped 
oscillations, see Figs. [5] and [31 The results shown in Figs 
[21 m and [5l are obtained by numerical evaluation of the 
sums ([T^l) and pU)) with the second form of the scalar 
product p3l) . For finite T/eq, the motion of the gas cen- 
ter of mass is not strictly a relaxation to the asymptotic 
value zg(oo). The distance zg(^) ~ zg{oo) reaches values 
of the order of Z for times t 3> {Z/gY^'^. However, these 



FIG. 4: Gas density as a function of height for fermions, ^ = 
0.1, Z = T/mg, T = 50eo and times t = 1.3{2Z/gY^'^ and t = 
2.1{2Z/ gY ■ These times correspond to the first minimum 
and the second maximum of the gas center of mass, see Fig. 
O The full lines are the large A^ density at the same times. 
The dashed line is the large A'' density at infinite time. 



fluctuations disappear in the limit T 3> co- For example, 
for C = 10 and Z — \, ZQ(t) — zg(oo) is of the order of 
0.04Z for t ~ 40(Z/.g)i/2 at T = 5eo but reaches again 
values of the order of 0.02Z only for t p; 110(Z/g)i/2 at 
T — lOeo- The dependence on C, and Z of the center of 
mass time evolution (|22p is not simple, see Fig(31 We ob- 
serve that the number of bounces increases with Z and 
that the boson and fermion curves are very close to each 
other already for C = 0.1. 

To obtain the large N limit of the density (|10p . we 
write it in a form similar to p^ - lj^O)) . In this case, the 
terms p = q are present in the sum. Their contribution to 
the analog of the function J is proportional to S{u!) and 
hence gives rise to a static component of the gas density 
wich we denote by poo ■ The other terms can be evaluated 
with the help of the approximations used above to derive 
We find 



cs-kX 



p{z,t) 



s>0 




dx I I 9 

—^ cos X\ - — t 



sin ^y(s^^Y — x"^ cos ■^/(sTr)^ — zx'^/Z 



(23) 



where X ~ min(l, ^ Z/ z). For the asymptotic particle 
number density, we obtain 



Poo[z) 



1 



V^A 



>l/2 



dx 0(x — mgz/T^ [x 

V^ ^x - mgz/T e^"^/C T 1 



.(24) 



The large N expressions (|23| and (|24ll agree very well 
with the exact results already for T/cq = 50, see FigHl 
For Z = 0, dMl) simplifies to ^ with Z = 0. In this 
case, the density p remains constant and equal to the 
thermal equilibrium density profile determined by the 
Hamiltonian Hq, energy E and number N. For Z ^ 0, 
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the populations |^aP remain constant. In this case, the 
effect of the bath is generically to destroy the coherences 
and the gas state turns asymptotically into 



^~ = IIl*'^l'l«)(°l 



(26) 



The corresponding particle number density Pqo{z) is 
equal to the constant component of ([4]). In the general 
case, [Hi,Ho] ^ 0, provided the coupling to the bath is 
weak enough, the asymptotic gas state is 



^c. 



= z- 



^-EJTb 



a){a\ 



(27) 



FIG. 5: Gas density as a function of height for bosons, con- 
densate fraction of 0.1, Z = T/mg, T = 50eo and times t = 
and f = 2(2^/3)^/^ 



the asymptotic density profile is very different from the 
thermal equilibrium one. For z 3> T/mg, it coincides 
with the initial density (fTT)) . For z < Z, the 8 function 
in ([M)) can be replaced by 1 and the derivative dzPoc 
is obviously positive and diverges for z — Z. The den- 
sity Poo thus assumes a maximum at a height z > Z, 
see Fig. |4l This conclusion also applies to the 3D case 
mentioned at the end of llllBI For a Bose gas below the 
condensation temperature, the behavior of the conden- 
sate is clearly different from that of the thermal cloud, 
see Fig. [5]where the condensate fraction is 0.1. The par- 
ticle number density of the condensate at times multiple 
of 2{2Z/gy/'^ is essentially identical to its initial density 
[l65| whereas the non-condensed gas relaxes to its asymp- 
totic profile. 



V. INFLUENCE OF AN ENVIRONMENT 

Here, we compare the result of the relaxation process 
in the presence and absence of a heat bath. If the gas is 
coupled to a heat bath, the complete system consisting 
of the gas and its environment is described by a product 
Hilbert space Hgas (X" Henv and by a Hamiltonian of the 
form H = Hq + Hb + Hj where Hb is the bath Hamil- 
tonian and Hj the interaction between the gas and the 
bath. In other words, the gas can exchange only energy 
with its environment. Owing to the coupling Hi, the 
state of the gas, initially given by ([3]), becomes mixed. 
To discuss the infiuence of the bath, we rewrite the initial 
state ft= |*)(^| as 



n^Y.'^,n\a){b\ , M/,= ^ 



*. 



(25) 



\a)eHE 



where we denote by Latin letters the TV-particle eigen- 
states of the Hamiltonian Hq. Two cases must be distin- 
guished. In the special case where Hq and Hj commute, 
the gas Hamiltonian is a constant of motion and hence 



where Z — J2a 6xp(— -Eq/Tb), Ea is the eigenenergy cor- 
responding to \a) and Tg is the bath temperature |ll.[l4|. 
The gas relaxes into thermal equilibrium at the temper- 
ature of the bath. 

If the gas is isolated, its state remains pure but, as seen 
in the previous section, the particle number density can 
relax into the steady component p^o of ^. Therefore, 
the asymptotic density profile of the isolated gas is not 
modified by an energy-conserving environment, i.e, such 
that [Hj,Hq] = 0. However, the relaxation towards this 
density is obviously different and depends on the details 
of the Hamiltonians Hb and Hj. In the general case, 
if the coupling to the bath is weak enough, the density 
profile poa determined by the initial state |^) can be ob- 
served in a transient regime before the bath imposes its 
temperature on the gas. 



VI. CONCLUSION 

In this paper, we have studied an isolated perfect quan- 
tum gas confined by gravity and initially trapped above 
a certain height. The gas was assumed to be in a pure 
state of macroscopically well-defined energy and to evolve 
under Schrodinger dynamics. We have derived the ex- 
pression for the time-dependent gas density profile in the 
limit of a large particle number. This single-particle dis- 
tribution was found to depend not on the exact micro- 
scopic state of the gas but only on a few thermodynamic 
parameters: the characteristic height of the initial con- 
finement, the microcanonical temperature and chemical 
potential corresponding to the gas energy and particle 
number. Damped oscillations of the density profile were 
obtained. This relaxation behavior is a manifestation of 
the many-body nature of the system. The time evolu- 
tion of a single confined particle is radically different. As 
an important consequence of this difference, the time- 
dependent density of a Bose gas below the condensation 
temperature consists of two clearly distinguishable com- 
ponents corresponding to the condensate and the thermal 
cloud. However, though the gas density relaxes, in gen- 
eral, towards a well-determined profile, this asymptotic 
density is very different from the thermal equilibrium one 



determined by the gas energy and particle number. Con- 
trary to the latter, the former depends also on the initial 
height of the gas and presents a maximum. 

We have also derived general results about the time- 
evolution of an isolated many-body system. We have 
shown that, the single-particle distributions of a typi- 
cal pure state of macroscopically well-defined energy are 
identical to that of the microcanonical mixed state at 
the same energy. We stress that this equivalence holds 
for any many-body system, our proof does not depend 
on the absence or presence of particle-particle interac- 
tions. Moreover, as the single-particle distributions are 
the same for almost all pure states of a given macro- 
scopic energy, the microcanonical results should be valid 
for a large class of statistical ensembles and apply to ex- 
perimental measures obtained with different systems pre- 
pared " under identical experimental conditions" . A cou- 
pling of the system to environmental degrees of freedom 
modifies the time evolution of the single-particle distri- 
butions. However, as discussed in the previous section, 
the asymptotic distributions found assuming the system 
is truly isolated are not changed by an energy-conserving 
environment and could be observed in a transient regime 
if energy is exchanged with the environment. 

We finally comment about the particle-particle inter- 
actions. We first emphasize that our results show that 
these interactions are not necessary to obtain a relax- 
ation behavior of physically relevant degrees of freedom 
of a many-body system. As we have seen, this behavior 
results from the properties of the microcanonical single- 
particle density matrix. Furthemore, the dilute regime is 
experimentally accessible. In this regime, for Bose gases, 
the interactions cannot be neglected below the conden- 
sation temperature as the condensate size is essentially 
fixed by the confining potential but they might play a 
minor role above this temperature [2l|. The theoreti- 
cal study of the influence of interactions between the gas 
particles and other issues are left for future work. 
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where X = Z/zq. The first equality is obtained with the 
help of the equation [d^ — x]A{x) = and the second 
one using this equation and integrations by parts. As 
A{aq) — 0, the above result for X — gives the orthogo- 
nality of the wavefunctions (f>k and (jfq and, for k — q, the 
equality of the derivatives with respect to A at A = of 
both sides gives the normalization constant of (f>k • 

By integrations by parts and using the fact that (j)q 
are the one-particle eigenfunctions of the Hamiltonian 
([1]) with Z = 0, the time derivative of the gas center of 
mass (TTSl) can be written as 



dfZG 



I 
Nm 



E' 



CpCq)^ 



it{e 



dz(j)q{z)dz4>p. (30) 



This expression is valid for any confining potential and 
simply states that the velocity of the gas center of mass 
is the average momentum of the gas divided by its mass 
Nm. In the same manner, we find 



d' 



t^G^T^ 



9 x^> 






it(ep — Cg) 



(31) 



where we have used the normalization of the density (jTC 
and 9z(/)p(0+) = (2m^g)^/^ for any p. As {c\,c^ = {c^qC^ 
is real, Y.p^qk(%c.<^ I k^v " '^9) = and dtZaiW) = 0. Con- 
sequently, the center of mass zq is given by (|19p . 

At initial time, the gas density vanishes at z = and 
hence we expect an acceleration 5^zg(0) = —g. This 
can be shown in the following way. Let us consider a 
wavefunction ^(^z) which vanishes for z < h where h > 0. 
We expand it on the basis {(f>p} and write its derivative 
at z = as 



To evaluate the scalar product P^ . we write 

dxA{x + aq)[d'^~x]A{x-X + ak) (28) 



X 



(flfc - A) / dxA{x + aq)A {x ~ X + a^) 



X 



^A{X + aq)d^A{ak) 



+ aq dxA{x + aq)A{x — X + Uk) 

Jx 



(29) 



dMO) = (2^^5)1/2 / dzifiiz) V 4>p{z) = 0. (32) 



We deduce from this property of the sum ^ 4ip{z) that 

^p,q(^l^q) = where {clcg) is given by ^^ and ^^. 

Therefore, J2p^q{4^q) = -J2p{4^p) = ~N and we ob- 
tain the expected initial acceleration. 
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